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Abstract —This paper studies an energy efficient design of pre¬ 
coders for point-to-point mnltiple-input-multiple-ontput (MIMO) 
orthogonal frequency-division multiplexing (OFDM) systems. 
Differently from traditional approaches, the optimal power al¬ 
location strategy is studied by modelling the circuit power as 
a rate-dependent function. We show that if the circuit power 
is a constant pins an increasing and convex function of the 
transmission rate, the problem of minimizing the consnmed 
energy per bit received can be reformulated as a convex fractional 
program and solved by means of a bisection algorithm. The 
impact of the some system parameters is investigated either 
analytically or by means of computational results. 

I. Introduction 

Nowadays, energy efficiency for wireless communications 
is becoming a main economical and societal challenge IT]. 
Hence, instead of maximizing the information rate under a 
certain power constraint, engineers wish to maximize the 
transmission rate per consumed power. This is also referred to 
as the energy efficiency (EE) maximization problem. Among 
the early contributions in this area, the authors in |2| originally 
considered the transmission rate per consumed energy to 
adapt the transmission mode (i.e. number of streams, space 
time signaling, MIMO detection, etc.) of a MIMO-OEDM 
system. Paper 0 investigates the optimal precoding strat¬ 
egy for EE maximization in MIMO systems, showing that 
without considering the power used to feed circuitries, the 
optimum transmission power will be zero for a Gaussian 
MIMO channel. In ID, the authors used the mathematical 
tool of convex fractional programming to address the power 
allocation problem for EE maximization in OEDM systems, 
showing the optimality of the waterfilling power allocation 
policy. More recently, the authors of 0 extended this work 
to a multi-user and multi-cell scenario. 

An alternative formulation of the original EE maximization 
problem was proposed by the authors of 0 who stated it as the 
minimization of the consumed energy per bit of information 
received. This route has also been pursued in Q and 0 to 
investigate the power allocation problem in OEDM systems, 
confirming the waterfilling nature of the optimal solution. 
All the papers above, however, consider the overall circuit 
power consumption (except the transmit power) to be either 
zero or a constant value. On the contrary, the authors of m 
accounted for the variability of the circuit power consumption 
for discrete-time base-band signal processing (due to different 


modulation orders and coding rates) by modelling this power 
as being dependent on the transmission rate. 

In the current work, we consider the EE design of precoders 
for MIMO-OEDM systems, formulated as the minimization 
of the consumed energy per bit of information received. As 
a main result, capitalizing on the quasi-convexity of the EE 
function, we derive the optimum precoding and power allo¬ 
cation strategy. Eventually, the impact of system parameters 
is mathematically analyzed and the results are confirmed by 
means of computational results. 

In this paper, boldface lowercase letters and boldface up¬ 
percase letters represent vectors and matrices, respectively. 
Notation {■)^ and tr (•) denote conjugate transpose operation 
and trace of a matrix respectively. | • | is the determinant of 
the matrix. (A)^ denotes the (m, m)-th entry of the matrix A. 
In represents the NxN identity matrix. (j>'{x) stands for the 
first-order derivative of the function (j){x). Xds is X expressed 
in decibels. The notation [xjjj" denotes max {0, x}. Einally, we 
denote the expectation operation by E{-}. 

II. SYSTEM MODEL 

Consider a MIMO-OEDM system with M transmit anten¬ 
nas, N receive antennas and K subcarriers. Assuming perfect 
synchronization, the observation model at the fc-th subcarrier 
is given by 

yfc = V^HfeXfe-fiifc. (1) 

Hfe S i^NxM j-jjg MIMO channel matrix for the /c-th 
subcarrier, whose entries are assumed to be independent and 
identically distributed (i.i.d.) zero-mean circularly-symmetric 
complex Gaussian random variables with unit variance. /3 = 
is the path loss budget where n is the path loss 
exponent, Gi is the gain factor at d = Im and Mi is the 
link margin accounting for the hardware process variations and 
other noise and interference ifTOI . Xfe is an M x 1 vector which 
denotes the transmitted symbols with a transmission power 
Pk = E{x{{xfc}. Hfc is a zero-mean complex additive white 
Gaussian noise vector with covariance E{nfcn{{} = (t^In, 
where tr^ = BNqN f, B is the bandwidth per subcarrier, Nq 
is the one-sided noise power spectral density, and Nf is the 
noise figure defined as in Qo]. 

According to i), the total power consumption of the 


transceiver can be modeled as 
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where lo is the efficiency of the power amplifier, P^, = 
PtcM + PrcN accounts for the power needed to feed the 
radio frequency (RF) chain ifTOl . while K(j) J2k=i is the 
rate-dependent term due to the discrete-time basebancl signal 
processing, k represents a constant coefficient, 9k is the infor¬ 
mation rate at the fc-th subcarrier, and (/){■) is assumed to be 
monotonically increasing function with (j){0) = 0. Tradition¬ 
ally, the EE function is defined as the information rate divided 
by the total energy consumption: 0^^ /-Ptotai- in 

this paper, we pursue a different route by investigating the 
consumed energy per bit of information received, i.e. 


e = 


hJ:tiPk + Pc + ^^^ {pEtiOk) 


{Joule/bit). (3) 


In the following, we will see that thanks to this formulation 
of the performance metric, we can restate the problem as a 
convex program and derive the optimal solution. 


III. Energy efficient precoder design 

In this section, we first restate the original fractional pro¬ 
gram as a convex optimization problem and then we derive the 
optimal precoding strategy. Einally we propose an algorithm 
achieving the minimum energy consumption per bit received. 


A. Precoder optimization 

Our aim is to minimize the function defined as in ([2i. The 
problem can therefore be formulated as 


min 

{QO 
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where Qj, is the covariance transmission matrix for the fc-th 
subcarrier with Pk (Qfc) = tr (Q^) = Ejx^Xfc}. Assuming 
Gaussian codewords, we have 
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We first show the following lemma. 

Lemma 1: The optimal structure of Qk is = V^Q^V^, 
k = where Q;, is a diagonal matrix and Vfc comes 

from the singular value decomposition (S VD) of Hfc that = 

UfcSfcVf. 

Proof: Assume that is the solution of @ and 
V^Q^Vfc is not diagonal. According to SVD of H^, we have 

/3Hf Hfc = VfeAfeVf, (6) 


where Ufc is a unitary matrix and = P'Ejf'Ek is a diagonal 
matrix. Define Q^, = V^Q^Vfe, where tr = tr (Q^), 


thanks to Hadamard’s inequality, we have 
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where the first equality is due to Sylvester’s determinant 
theorem that log 2 |Ip + AB I = log2 |Ig + BA|, where A is 
a p X q matrix and B is a g x p matrix. Therefore there exists 
a diagonal matrix Q;, such that 


Pk = tr < tr = Pf, 


(8) 


and 
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Then from Q, defining Q'^, = U^Q^lIfc, we have 

e{Qt,...,Q'fc,...,Qj^} <e{Qt,...,Q^...,Q^^}, (10) 

which is contradictory with the fact that minimizes the 
objective function. Therefore, every matrix Q^., k = 1,K, 
should be a diagonal matrix. ■ 

In the following, we only consider the subset of matrices Qj. 
that have the structure as in Lemma 1. Hence, the rate for 
subcarrier fe can be rewritten as 
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which defines 6k,m- As 


and 


it appears that there is a one-to-one mapping between the rate 
and the covariance matrix, and that Pk can be written as a 
function of 9k,i, ■■■, 9k, m- Therefore, instead of formulating the 
problem in the variables Q;,, we can restate it in the variables 
9k,Tn, which IS similar to the approach of 0: 
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where, for a given t, this problem is convex. The corresponding 
Lagrangian is: 
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which, after some manipulation of (fTSI l. leads to 
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Substituting ( fT9] l into (fT2l i get 
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Proof: Let us denote p(0) = + K(j){BQ). We have 

Pc + 5 ( 0 ) 


e(0) = 


B 0 


(24) 


(15) 


Therefore, the optimal solution must fulfill the following KKT 
conditions: 
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Defining the water level as being fj, = the next step is 

to find its optimal value. This can be obtained by rewriting 
the total power (P = ^k=i ^k) as a function of the total rate 
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where L is the number of eigenchannels receiving a non¬ 
zero power among the min{M, TV} x K space/frequency 
subchannels. pi denotes the transmission power and A; denotes 
the channel gain for the Lth subchannel. The derivative of the 
transmit power with respect to the total rate is given by 


Assume that 0 1 < 02 B^nd 0 ^ A <C 1. If 

e(A0i + (l-A)02) >e(0i), (25) 

we have 

Pc + g (A 01 + (1 — A) 02 ) Pc + g(0i) 

B(A0i + (1-A)02) BQi ' ^ ’ 

Noting that, thanks to the convexity of p(0), 

g (A 01 + (1 - A) 02 ) < Xg (0i) + (1 - A) 5 ( 02 ), (27) 

we get 

PcB (A01 + (1 - A)02) + Bg{&i) (A01 + (1 - A)02) 

< PcBQi + BQig (A01 + (1 - A)02) 

< PcBQi + BQiXg (0i) -f B0i (1 - A) <? (02) (28) 

which after some manipulations, leads to 

ft < (29) 

02 — 01 

If instead of (l25] t, we assume that 

£ (A01 + (1 — A)02) > ^(02), (30) 

similarly we get 

„ ^ 013 ( 02 ) - 023 ( 0 i) 

02-01 ■ 

Equations ( l29b and OTl i cannot be met simultaneously, which 
means that 

£ (A 01 + (1 — A) 02 ) < max {e( 0 i),e( 02 )} , (32) 

which implies that e( 0 ) is a quasi-convex function. ■ 

Lemma 2 means that + K(j){BQ) should be convex to 
guarantee that the objective function is quasi-convex. 

The optimal transmission rate should fulfill the following 
condition: 

£'(©) = ^ (^)'otal(0) - = 0- (33) 

Let us consider a generic convex model for the baseband 
power consumption, in the form of <j){BQ) = {BQ)°‘ (a > 1). 
Then, we obtain 
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from which it can clearly be seen that P(0) is a strictly 
increasing and strictly convex function. 

Then we have the following lemma: 

Lemma 2: If +K(f){BQ) is convex, then e(0) is quasi- 


where we made explicit the dependence of fi with respect to 
0. Hence, the value of the total information rate minimizing 
the consumed energy per bit must be such that 
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B. Algorithm description 

In the previous subsection, we showed that the original frac¬ 
tional program can be solved by simply finding an information 
rate value that satisfies dlSl) . Unfortunately, such a value cannot 
be obtained in a closed form. Hence, we resort to an iterative 
approach based on a bisection algorithm. The algorithm is 
provided hereafter. The initial value 0° can be set to any 
arbitrary strictly positive value. The iteration stops when the 
search interval becomes smaller than a prescribed threshold 5. 

1: Set 0; = 0; Ou = 0° 

2: Calculate /lu and Pu at the point 0^. Calculate e'(0u). 

3: while s'{Qu) < 0 do 
4: Bu = Bu X 2 

5: end while 

6: while Bu — 0; > <5 do 
7: 0c = 0.5(0;-h 0u) 

8: Calculate n(0c) 

9: if e'(0c) = 0 then 

10: go to line 15 

11: else if n(0c) > 0 then 

12 : 0 „ = 0 c 

13: else 

14: 0; = 0c 

15: end if 

16: end while 

TABLE I: Algorithm description 

IV. Impact of system parameters 

In this section, we analyze the impact of some system pa¬ 
rameters on the EE link performance metric. We first consider 
the impact of the distance d on the minimum consumed energy 
per received bit when all the other parameters are kept to a 
given constant value. 

Assume two different values of the transmitter-receiver 
separation d, namely di and c? 2 , where d2 > di. Denote by 
e(0; d) the consumed energy per bit evaluated for a rate 0 at 
distance d. Denote by Pc{d) the power consumption needed 
to feed the RE chains for distance d and by 

0*(d) = argrnin{e(0; d)} (36) 

the value of the rate 0 corresponding to the minimum con¬ 
sumed energy per bit at distance d. Einally, denote by P(0; d) 
the total transmit power evaluated at distance d to achieve the 
rate 0. Then we have 

e(0*(d2);d2) 

+ Pc(d 2 ) + iBe*{d2)) 

B&*{d2) 

+ Pcidi) + ncjr (P0*(d2)) 

^ P0*(d2) 

^ +p^(d^) + «<^(B0*(d2)) 

^ P0*(d2) 

= e(0*(d2);di) > £(0*(di);di)- (37) 

The second equality comes from (|6]), (l2Tl i and (l22T i observ¬ 
ing that, when the distance d increases from di to ^ 2 , the 


channel gain of every subchannel (A; in (ItTIi ') is multiplied 
by (d 2 /di)“". Hence, according to (l22l) . the rate can be kept 
constant if the water level is multiplied by (d 2 /di)". Since the 
rate is a strictly monotonic function of the water level, this is 
the unique possible choice. Thus according to the waterfilling 
policy, the number of subchannels with non-zero power is 
not changed, and the total transmit power is multiplied by 
((i 2 /di)". Pc is independent of the value d. The last inequality 
is due to the fact that 0 *(di) minimizes the function when 
d = di- Therefore we have, £( 0 *(d 2 );d 2 ) > e(0*(di);di), 
meaning that the minimum number of Joules per received bit 
increases with the distance. 

Let us now analyze the impact of distance d on the rate value 
at the optimal point. Let us define Ptotai(0:'7) as the total 
power written as a function of 0 and d. With the definitions 
provided above, we have the following theorem; 

Theorem 3: If </)(•) is linear, then 0 *(d 2 ) < 0*(di) for 
d 2 > di. 

Proof: Since 0*(di) is the optimal point that minimizes 
the objective function when d = di, according to (l33T l. 


e'( 0 *(di),di) = 


P' (0*(dl) d,) -Pt°tai(e*(di),di 
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If </>(•) is linear, after some manipulations, we get 
P'(0*(di),di) :M^iMl)+P,(di) 
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where the derivative of the power in the first line together 
with equation (|2^ yields to the third line. Using the fact 
that e( 0 ,d) is quasi-convex, e'( 0 ,d) > 0 for 0 > 0 *(d), 
£:'( 0 ,d) < 0 for 0 < 0 *(d), which is proven in detail in 
Theorem 1 in m. Therefore we have 0*(di) > 0 *(d 2 ). ■ 

A similar analysis can be carried out for parameters cr^, 
Ptc and prc- We however omit it for the sake of concision. It 
appears that increasing the noise variance has an impact 
similar to increasing the distance, meaning that the optimal 
rate value decreases with increasing noise variance. As far as 
Ptc and Pec are concerned, increasing their value also leads to 
a decrease of the corresponding optimal value of the rate. 
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Fig. 1: Power consumption and EE comparison for different 
d with M = N = 4 and K = 64. 
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Eig. 3; Power consumption and EE comparison for different 
antenna configurations with d=10m and K = 32. 
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Eig. 2: Power consumption and EE comparison for different 
K with M = N = A and d=50m. 


V. Numerical results 

In this section, we illustrate our analytical findings by means 
of numerical results. The system parameters are set as follows: 
the bandwidth for each subcarrier is set to B = lOKHz, 
Ptc = 82.5mW, Prc = 105.5mW HO), and we fix k = 5 x 10“® 
as in a. We also select the following values: n = 3.5, 
GodB = “ {GidB + ^idB) = —70dB where GidB = 30dB 
is the gain factor at d = Im and MidB = 40dB ifTOl . The 
noise power spectral density is set to N^dB = —170dBm/Hz, 
the noise figure to = lOdB/Hz and the amplifier 

efficiency is chosen to be w = 0.4. Unless otherwise specified, 
we use a parameter a = 1 for the rate dependent power 
consumption term. Einally, we set 5 = 0.01 as the tolerance 
of the bisection algorithm. The results are averaged over 1000 
different channel realizations. 


A. Effect of link distance d 

FiglU reports both the different power consumption terms 
and the value of the objective function e at the optimal point, 
for different values of d = 10,50,100m. The number of 
antennas and the number of carriers are set to M = = 4 

and K = 64, respectively. As said before, the value of the 
objective function e increases with the distance. So does the 
transmit power. As expected from Theorem 0 the optimum 
rate value decreases with increasing distance which means 
that the rate dependent power consumption also decreases with 
increasing distance. In view of these different dependencies of 
power consumption terms with respect to distance it is however 
difficult to predict the evolution of the total power with the 
distance. 

B. Effect of bandwidth enlargement 

Eig|2] reports both the different power consumption terms 
and the value of the objective function e at the optimal 
point, for a number of subcarriers K growing from 1 to 256, 
meaning different bandwidth sizes KB. We select d = 50m 
and M = = 4. As it is expected, it can be observed that 

the consumed energy per bit decreases when the number of 
subcarriers grows. As a matter of fact, when the number of 
carriers is increased from Ki to K 2 the solution space for 
Ki is contained in the solution space of K 2 and the solution 
for K = K 2 cannot be worse than the solution for K = Ki. 
It is also observed that both the transmit power and the rate 
dependent power grow with K, meaning that the rate increases 
with an increasing value of K. 

C. Effect of antenna configuration 

Einally, we investigate the impact of antenna configurations 
on our objective function. Eigs. [3 and |4] report the value of the 
objective function at the optimum for 3 different sets of config¬ 
urations. In the first set M and N take simultaneously values 
of {2,4,6, 8}, which is beneficial for the spatial dimensions 
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Fig. 4: Power consumption and EE comparison for different 
antenna configurations with d=50m and K — 32. 
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Fig. 5; Power consumption and EE comparison for different 
antenna configurations with d = 10m and K = 32. cj){BQ) = 
which is nonlinear. 


captured. In the second set M = 2 and N G {2,4, 6 , 8}, which 
corresponds to receive diversity. The third set corresponds to 
transmit diversity; M G {2,4, 6 , 8 } and N = 2. As it appears 
from figure [2 for the system parameters selected here and 
with d = 10m and K = 32, set 1 is more beneficial when the 
number of antennas increases. This is mainly due to the pre¬ 
log factor associated with the spatial multiplexing gain. On the 
contrary, it is detrimental to only increase M (the number of 
transmit antennas) or N (the number of receive antennas). This 
can be explained by the fact that the short distance between 
transmitter and receiver translates into a small transmit power, 
which turns out not to be the dominating term in the total 
power consumption. Otherwise stated, the total power is more 
impacted by the constant term, meaning that the optimal point 
corresponds to a high SNR value 0 . Hence, when the number 


of transmit (receive) antennas increases, the diversity gain 
cannot compensate the increase of power consumption due 
to the additional RE chains. Eig. 0] compares the three sets 
of antenna configurations for a larger distance, i.e. d = 50m. 
While the conclusion remains the same for the first set of 
configurations, increasing the number of antenna turns out 
to be beneficial for the two other sets, exploiting diversity. 
As a matter of fact, the transmitted power increases due to 
distance and becomes dominating in the total power. The 
optimal operating point corresponds to a lower SNR 0 at 
which diversity gain prevails over the additional power due to 
the increasing number of RE chains. 

Interestingly, our previous conclusions are sensitive to the 
value of a. Eig|5] reports result for a scenario where the 
(/>{■) function is nonlinear with a = 1.2. The first set of 
configurations, i.e. M = N, is considered. It now turns out 
that increasing M = N is detrimental, in opposition to what 
we had for a = \. This is due to the rate dependent power 
consumption that grows with the multiplexing gain faster than 
the information rate. This shows the crucial role played by the 
power consumption model in optimally designing the link. 

VI. Conclusion 

In this paper, we studied the energy efficient design of pre¬ 
coders for point-to-point MIMO OEDM systems. We showed 
that for the total power made of a constant term plus another 
one that is increasing and convex with the transmission rate, 
the consumed energy per bit is a quasi-convex function of 
the total transmission rate. Thanks to that, the problem of 
minimizing the consumed energy per bit could be reformulated 
as a convex fractional program and solved by means of a 
simple bisection algorithm. The effects of various system 
parameters on the optimal value of the objective function 
have been analyzed and illustrated by means of computational 
results. 
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